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SUMMARY 

This report describes an investigation to predict first-ply failure and final fracture in selected composite laminates 
subjected to inplane loads. The laminates were composed of glass fiber and graphite fibers in epoxy matrices. Fail- 
ure envelopes based on first-ply failure and laminate fracture were generated for combined loading of these lami- 
nates. Predictions were evaluated by micromechanics-based theory and progressive fracture. The results show that, 
for most cases, combined tensile loading significantly enhanced the laminate fracture stress in comparison to the 
uniaxial loading. 


INTRODUCTION 

One of the persistent difficulties in composite design and analysis has been the inability to predict laminate frac- 
ture under uniaxial and/or combined loading by using either unidirectional composite (ply) data or micromechanics 
with pristine constituent material properties. The difficulty has been compounded many times over by the existence 
of many and diverse failure criteria. Thus it became apparent that a formalized comparison of the various failure 
theories to the measured data would be instructive and very useful. 

To that end, A.J. Hinton (Head of the Advanced Composite Structures Team, DRA, Fort Halstead, UK) and 
P.D. Soden (Senior Lecturer in Mechanical Engineering, UMIST, UK) organized an exercise to confirm state-of- 
the-art methods for predicting failure in composites. World-renown researchers in the failure of composites were 
invited to participate. The organizers provided them with (1) a list of laminates, (2) unidirectional composite data, 
(3) pristine constituent material properties, (4) unidirectional composite stress-strain curves, and (5) recommended 
standardized formats to present their results. Results shown in the various figures in this report comply with the 
recommended formats. The objective of this report is to describe the results we obtained during our participation in 
that exercise. 

Our contribution to this exercise consisted of two major parts: (1) using the micromechanics-based uniaxial com- 
posite strengths described herein to generate first-ply failure envelopes for combined loading, and (2) using progres- 
sive fracture to generate laminate fracture envelopes and stress-strain diagrams for uniaxial and combined loading. 
The effects of residual stresses on laminate first-ply failure were evaluated, as were the differences between pristine 
and in situ constituent properties. The pristine constituent material properties were calibrated to the unidirectional 
composite properties so that the micromechanics used would reproduce the unidirectional composite properties 
provided. Details are described in this report, along with the specific theories and computer codes used. 

The laminate failure envelopes were generated by two computer codes: (1) the Integrated Composite Analyzer 
(ICAN) and (2) the Composite Durability Structural Analyzer (CODSTRAN). For completeness, these computer 
codes are briefly described in the following sections. 



ICAN COMPUTER CODE 


The authors have conducted extensive research at NASA Lewis Research Center, developing composite mechan- 
ics theories and analysis methods ranging from micromechanics to finite elements. These theories and analysis 
methods account for environmental effects and are applicable to intraply hybrid composites, interply hybrid compos- 
ites, and combinations therefrom. Most of these theories are represented by simplified equations and have been cor- 
roborated by experimental results and finite-element analysis. The composite mechanics theories, with their 
simplified equations, constitute a structured multiscale formalism that is (1) “upward integrated” (synthesized) from 
material behavior space to structural analysis and (2) “top-down traced” (decomposed) from structural response to 
material behavior space (see fig. 1). 

This theory has been incorporated into a computer code called ICAN (Integrated Composite Analyzer). The fol- 
lowing paragraphs give a brief history of the developments in composite mechanics and the related computer codes 
that led to the evolution of ICAN. 

The need for multilevel (multiscale) analysis in designing structural components made of multilayered fiber com- 
posites was recognized about 30 years ago (ref. 1). Shortly thereafter, a multilevel analysis computer code (Multilay- 
ered Fiber Composite Analysis (MFCA)) was developed; it was found to be efficient in predicting the structural 
response of multilayered fiber composites, given the constituent materials properties, the fabrication process, and the 
composite geometry (see ref. 2). 

Since intraply hybrid composites were the logical sequel to conventional and interply hybrid composites, theoret- 
ical and experimental investigations were subsequently conducted on the mechanical behavior of intraply hybrids 
(refs. 3 to 5). The theoretical methods and equations described in references 3 to 5, together with those for hygro- 
thermal effects (ref. 6), have been integrated into a computer code that predicts hygral, thermal, and mechanical 
properties of intraply hybrid composites, thereby enabling their “design.” The computer code is called “INtraply 
HYbrid Composite Design,” or INHYD (ref. 7). 

The ICAN computer code is a synergistic combination of MFCA and INHYD, along with several significant en- 
hancements. It utilizes the micromechanics embodied in INHYD and the laminate theory of MFCA to build a com- 
prehensive analysis and design capability for structural composites. ICAN also offers the following unique features: 

(1) Ply stress-strain influence coefficients 

(2) Microstresses and microstress influence coefficients 

(3) Stress concentration factors at a circular hole 

(4) Predictions of probable delamination locations around a circular hole 

(5) Poisson’s ratio mismatch details near a straight free edge 

(6) Free-edge stresses 

(7) Material cards to be implemented as input for finite-element analysis by several commercial computer codes 

(8) Laminate failure stresses based on first-ply failure and fiber-breakage criteria, with and without hygro- 
thermal degradation 

(9) Transverse shear stresses and normal stresses 

(10) Explicit specification of interply layers 

In addition, ICAN has its own (dedicated) data bank of material properties for commonly used fibers and matri- 
ces, which can be easily modified by the user to add new constituent materials. 

ICAN is primarily designed to analyze the hygrothermomechanical response and the properties of fiber- 
reinforced/resin-matrix layered composites, given the local membrane loads and bending moments. Three types of 
layers are recognized by the program: (1) a layer in the standard composite system, which consists entirely of a pri- 
mary composite made of one type of fiber and matrix; (2) a layer in the intraply hybrid composite system, which 
consists of a primary composite and a secondary composite arranged in a prescribed manner within a layer (for pur- 
poses of identification, the primary composite in the hybrid is the one that constitutes the largest volume ratio); and 
(3) the interply layer, which consists of a matrix only. 

Complete details of the equations in the code are given in reference 8. The composite hygrothermomechanical 
properties can be predicted by using the micromechanics theories mentioned previously. Laminate properties are 
obtained with macromechanics and laminate theory, and local stresses and strains are determined from classical 
laminate theory. The free-edge stress calculations are based on the approach outlined in reference 9, with enhance- 
ments to accommodate the interply layer and the local characteristics of the adjacent ply. The stress concentration 
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factors around a circular hole are obtained by using the equations given in reference 10. Laminate failure stress 
analysis utilizes two criteria: (1) first-ply failure based on maximum strength and (2) first-ply failure based on fiber 
breakage. Complete laminate failure analysis uses two different ply combined-stress failure criteria and one interply 
delamination criterion for each specified load condition. 


CODSTRAN COMPUTER CODE 

CODSTRAN is an integrated computer code based on a new approach to progressive fracture of composite lami- 
nates and structures. It evolved from the following observations and/or concepts: (1) Any laminate or structural 
component can sustain a certain amount of damage prior to structural fracture (collapse). (2) During damage propa- 
gation, the laminate or component exhibits progressive degradation of structural integrity as measured by global 
structural behavior variables, such as loss of frequency, loss of buckling resistance, or excessive displacements. 

(3) The critical damage can be defined as the amount of damage beyond which degradation of the laminate struc- 
tural integrity is very rapid, having been induced by either a small amount of additional damage or a slight increase 
in loading. (4) Laminate damage is characterized by the following sequence: damage initiation, growth, accumula- 
tion, stable or slow propagation (up to the critical amount), and unstable or very rapid propagation (beyond the criti- 
cal amount) to collapse. 

These concepts are fundamental to developing formal procedures to identify the different stages of damage, quan- 
tify the amount of damage at each stage, and relate the degradation of global structural behavior to the amount of 
damage at each stage. 

CODSTRAN is a combination of composite mechanics (ICAN, ref. 8) and the finite-element method MHOST 
(ref. 1 1). It permits a simulation of local behavior as well as global structural behavior through an integrated com- 
puter code, which is shown schematically in figure 1 . The lower part of this figure describes the constituent proper- 
ties as functions of environmental and mechanical loading conditions. The criteria can be implemented by knowing 
the developed ply stresses, and ply strength can predict the damage initiation, growth, accumulation, and propaga- 
tion. As a result, the constituent properties are updated at every load increment. For example, if the ply longitudinal 
stress exceeds the allowable strength, then the fiber/matrix longitudinal moduli are replaced by negligible values so 
that, in essence, the ply does not carry any load and the stresses are redistributed to the surrounding plies. However, 
if a particular ply’s transverse strength exceeds its allowable strength, then only the matrix is assumed to have failed 
and the matrix modulus is replaced with a negligible value. Once the current constituent properties are determined 
(as in the left part of fig. 1 ), repeated applications of micromechanics, macromechanics, and laminate theory are 
used to assemble the global structural stiffness matrix, which is fed into the finite-element analysis. Thus, the left 
part of the figure depicts the integration (synthesis) of local damage conditions to global structural behavior 
(response). 

The nodal stress resultants are the outcome of the finite-element analysis. They are used to decompose the 
changes in the global response (e.g., laminate stresses and strains that result from incremental load or stiffness 
changes) on the local (micro) material stress/resistance. The load is incrementally increased only if there is no fur- 
ther damage due to changes in ply level stresses. Otherwise, only the material properties at the constituent level are 
updated at every iteration until a balance is reached between the applied loading and the local response. Overall 
structural equilibrium is maintained by iterations around the “cartwheel” until a specified convergence is reached. 
This procedure is illustrated in figure 2. The final result in terms of load as a function of global displacement is 
shown in figure 3. The schematics in figures 1 to 3, collectively summarize the fundamentals and implementation of 
this computational simulation of structural fracture in composites and composite laminates. 


MODIFICATION OF CONSTITUENT PROPERTIES 

In order to take full advantage of the simulation capabilities of ICAN and CODSTRAN, and since both of these 
computer codes are micromechanics (constituent materials) based, we modified the constituent materials provided 
by the organizers (tables I and II) so as to match the predicted composite unidirectional (ply) properties as shown in 
table III. A good reason for basing the ICAN/CODSTRAN simulations on micromechanics is that the in situ proper- 
ties of the constituents generally differ from those of their pristine state (ref. 1 2). The modified constituent properties 
we obtained and used are shown in tables 1 and II. A comparison of values for the matrix properties (see table I) 
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shows a substantial difference in matrix strength properties. Overall, the tensile strengths are about 50 percent lower, 
whereas the compressive strengths for the DFVLR matrix are about 30 percent higher. The shear strengths are also 
higher. The interesting point is that if the in situ matrix retains its pristine tensile strength, then the substantial tensile 
strength reduction could be caused by cracking due to thermal mismatch during the processing. The fiber properties 
(see table II) are about the same for both the as-provided and the calibrated cases. The unidirectional composite 
properties (table III) are almost identical. 


ICAN SIMULATION RESULTS 

Herein we summarize the simulation results obtained for biaxial stress failure envelopes by using the ICAN- 
modified constituent properties (table I) and first-ply failure. 


Unidirectional Laminates 

Biaxial failure stress envelopes were determined for a number of laminates. Figure 4 shows the envelope of an 
E-Glass/LY556 unidirectional composite subjected to transverse and shear loading (a y versus y). Since this failure 
envelope is symmetric about the cy axis, only the upper half is shown. The biaxial failure stress envelope for the 
T300/BSL914C unidirectional laminate subjected to longitudinal and shear loading ( g x versus x ) is shown in 
figure 5. It too is symmetric about the a^-axis, so only the upper half is plotted. 

The biaxial failure stress envelope for the E_Glass/MY750 unidirectional laminate subjected to biaxial longitudi- 
nal and transverse loading (g x versus G y ) is shown in figure 6. This failure envelope shows discontinuities from the 
tension-compression to tension-tension quadrants and from the compression-tension to compression-compression 
quadrants. The discontinuities are accentuated because of the large difference between the longitudinal and trans- 
verse tensile and compressive strengths. 


Angle-Ply Laminates 

The biaxial failure stress envelope for the E~Glass/LY556 [90°/±30°] s laminate subjected to a biaxial loading (G y 
versus g x ) stress state is shown in figure 7. Note that this failure envelope was generated by first-ply failure. This 
laminate is strong in compression, but relatively weak in tension. 

The biaxial failure stress envelope for E-Glass/LY556 [90°/±30°] s subjected to biaxial loading (x versus o x ) is 
shown in figure 8. This envelope, which is based on first-ply failure, is symmetric about the cyaxis, so only the 
upper half is plotted. The laminate is weak in both tension and shear, but strong in a compression. 

The biaxial failure stress envelope for the E-Glass/MY750 [±55°] s laminate subjected to G x and G y combined 
loading is shown in figure 9. This laminate is strong in ct^ tension in the cy g x compression-compression quadrant. 

The biaxial failure stress envelope for the AS4/epoxy 3501-6 [907±45°/0°] s laminate subjected to biaxial loading 
(G x versus G y ) is shown in figure 10. Note that this laminate is ultrastrong in the compression-compression quadrant, 
with the different types of combined loadings having a substantial synergistic effect in this quadrant. In the tension- 
tension quadrant the laminate is moderately strong, whereas in the second and fourth quadrants, it is weak. 


Effects of Thermal Stresses 

All the laminate failure envelopes were generated without accounting for residual stresses from lamination. These 
stresses result from the curing temperature differential and the thermal expansion coefficient (TEC) mismatch. It is 
instructive, therefore, to evaluate the effects of residual stresses on at least one failure envelope. Figure 1 1 shows the 
failure envelope for the AS4/epoxy 3501-6 [90°/±4570°] s replotted in full scale both with and without residual 
stresses. The residual stresses have a dramatic effect on the laminate biaxial failure stress envelope in the three quad- 
rants that have tensile stresses. This is an expected result because the plies in quasi-isotropic laminates have rela- 
tively high transverse ply stresses due to curing temperature differential and TEC mismatch. It should also be noted 
that resistance increased in the compression quadrant. 
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CODSTRAN SIMULATION RESULTS 


CODSTRAN simulation results are presented here for one unidirectional composite (for comparison purposes 
only) and for all the angle-ply laminates. 


Unidirectional Laminate 

The failure envelope generated by plotting CODSTRAN results is shown in figure 12. The plate was 10 mm wide 
by 10 mm long by 1 mm thick. The biaxial failure envelope was generated for versus O y . Although there are not 
enough points to obtain a smooth function, the points shown are exactly the same as in figure 4. This was an 
expected result since both CODSTRAN and ICAN use the same combined-stress failure criterion, which is stress 
based (ref. 13). From this we can conclude that linear composite mechanics is sufficient to develop combined-stress 
failure envelopes of unidirectional laminates subjected to combined loading. 


Angle-Ply Laminates 

Two combined-stress failure envelopes predicted by using the progressive fracture feature in CODSTRAN are 
shown in figure 13 for the E-Glass/LY556 [90°/±30°] s laminate subjected to combined a x versus ct , loading (as in 
fig. 7). The inner envelope is the initial damage failure envelope, and the outer one is the laminate fracture envelope. 
In figure 13 the initial damage envelope includes residual stresses, whereas in figure 7 it does not. Clearly, the 
residual stresses significantly affect the first-ply failure in quadrants with tensile stresses but do not affect the com- 
pression-compression quadrant, and the tensile stresses that cause laminate fracture are about 10 times those that 
cause first-ply failure. From this, two important conclusions may be drawn: (1) multi-ply laminates designed for 
first-ply failure under a x versus a y (tension-tension loading) will be significantly overdesigned; and (2) multi-ply 
laminates subjected to monotonic single loads to fracture must be designed by using progressive fracture analysis for 
cost-effective composite utilization. 

The CODSTRAN-simulated combined-stress failure envelopes for the Pi— Glass/LY 556 [90°/±30°] s laminate sub- 
jected to combined (x^ versus o x ) loading is shown in figure 14. The CODSTRAN simulations include residual 
stresses, which are significant contributors to first-ply failure. Comparing the inner failure envelope in figure 14 
with that in figure 8 shows that residual stresses decrease the laminate’s first-ply tensile and shear failures by a fac- 
tor of about 3. In addition, figure 14 shows that (1) progressive fracture (suitable for monotonic load) dramatically 
expands the tensile quadrant of the laminate failure envelope by a factor of about 15; and (2) the progressive fracture 
combined-stress failure envelope is not a smooth function because (a) the stress is redistributed following a local 
fracture, (b) the various failure modes are discrete and, therefore, not continuous, and (c) the jumps occur as the 
fracture progresses from one failure mode to the next. In summary, the important facts are (1) the laminate reserve 
strength from first-ply failure to laminate fracture is substantial when the laminate is loaded in combined shear and 
tension and (2) the laminate combined-stress fracture envelope is not a smooth function. 

Two laminate failure envelopes for combined o x and O y loading of the E-Glass/MY750 [±55°] s laminate are 
shown in figure 15. The inner envelope represents initial damage, and the outer one, laminate fracture. Interestingly, 
the two envelopes are practically identical when a x and a are both tensile. The observations and conclusions made 
in connection with figures 13 and 14 apply to figure 15 as well. 


Stress Versus Strain 


In this section we describe stress versus strain, as determined by progressive fracture via CODSTRAN. The 
stress-strain curves for all cases except one were obtained without including the effects of residual stresses. For that 
one case, the effects of residual stress were included. 

The stress-strain curve for the E-Glass/MY750 [±55°] s laminate is shown in figure 16. This laminate exhibits 
bilinear behavior: The first portion is linear up to a strain of about 0.5 percent; then the strain increases very rapidly 
up to 12 percent, or 24 times more, with a relatively small increase in stress (only 20 percent). The failure strain of 
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about 12-percent may be misleading because it occurs with relatively small load increments. In other words, lami- 
nate fracture was imminent past the 0.5 percent strain. One can infer from this that laminate configurations can exist 
with large strains before fracturing. Keep in mind, however, that this large strain is useful only for monotonic load- 
ing to fracture. The important conclusion is that progressive fracture is a suitable method for generating laminate 
stress-strain curves to fracture. 

Two stress-strain curves are shown in figure 17, where c y is plotted versus and e x for the E-Glass/MY750 
[±55°] s laminate: The strain is along the y-load direction, and the e x strain is along the x-load direction. Both 
curves are practically linear until fracture, with the z x strain exhibiting a minor deviation from linearity near the high 
strain region. It is interesting to note that the biaxial stress and strain reduced the fracture strain by a factor of 4 com- 
pared to the uniaxial loading shown in figure 16 while increasing the corresponding fracture stress by a factor of 6. 
Obviously, the stress-strain behavior of angle-ply laminates is dramatically different under combined loading in 
comparison to uniaxial loading. 

The stress-strain curve for uniaxial loading of the AS4/3501-6 [0°/±45 o /90°] s laminate is plotted in figure 18. 

This figure shows trilinear behavior. However, the last portion, beyond 1 -percent strain occurs with a very small 
increase in stress and happens very fast. The important feature is the nonlinear stress-strain behavior to fracture that 
the AS4/3501-6 [0°/±45°/90°] s laminate exhibits. The first linear portion ends at about 0.25-percent strain, which 
corresponds to trans-ply fracture in both the 90° and the ±45° plies. The second portion ends at about 1.0-percent 
strain, which corresponds to about 70 percent of the fracture strain of the fibers (see table II). 

The stress-strain curve for the AS4/3501-6 [0°/±45790°] s laminate subjected to combined-stress loading 
(a^/a^ = 2/1) is shown in figure 19. This stress-strain curve is linear to fracture at about 1 .5 percent strain, which 
corresponds to the fracture strain of the AS4 fibers (table II). The fracture stress is about twice that of the uniaxial 
case (fig. 19). This shows that some biaxial stress states substantially enhance the fracture stress of quasi-isotropic 
laminates. Unfortunately, the magnitude of the enhancement is load dependent and can be taken advantage of only 
for monotonic loadings to fracture. 

The stress-strain behavior of the E-Glass/MY750 [±45°] s laminate subjected to combined c } /a x = 1/1 loading is 
shown in figure 20. The useful fracture strain for this angle-ply laminate is about 2.5 percent, which is at the end of 
the first linear portion. Beyond this point the strain increases rapidly to about 6 percent with a relatively negligible 
increase in stress. The useful fracture strain of 2.5 percent corresponds to the strain just before the shear stress-strain 
curve becomes linear. Again we see that biaxial tensile loading enhances the fracture stress of [±45°] s angle-ply 
laminates, as it did for the [±55°] s angle -ply laminate (fig. 17). 

The stress-strain curves for the E-Glass/M Y750 [±45°] s angle-ply laminate subjected to biaxial a y f<3 x = 1/-1 
loading is shown in figure 21. Both stress-strain curves are linear to fracture at about 0.4-percent strain. This strain is 
considerably smaller than that in the tension loading case (fig. 20). The important conclusion is that fracture stress- 
strain of angle-ply laminates under biaxial loads with mixed signs (tension-compression) is substantially lower — 
only about 10 percent of that of the tension-tension cases. 

The stress-strain curve for the E-Glass/MY750 [0°/90°] s cross-ply laminate subjected to uniaxial loading 
(Xy/Oj = 0/1 is shown in figure 22. The stress-strain curve is linear to fracture at about 2-percent fracture strain, 
which is about 80 percent of the tensile fracture strain of the fiber (table II). And the fracture stress is about one-half 
of the fracture stress of the unidirectional laminate (table III). From this we may conclude that (1) cross-ply lami- 
nates loaded in uniaxial tension fracture at about one-half the developed tensile fracture stress of the uniaxial lami- 
nate and (2) cross-plying does not affect the linear stress-strain behavior of the unidirectional laminate. 


SUMMARY OF RESULTS 

An investigation was conducted to evaluate the first-ply failure and laminate fracture of selected laminates made 
from glass/epoxy and graphite/epoxy fiber polymer matrix composites that were subjected to uniaxial and combined 
loadings. The evaluation employed laminate theory for first-ply failure and progressive fracture for laminate frac- 
ture. Both are based on micromechanics and require constituent material (fiber and matrix) properties as inputs, as 
well as fiber volume ratio, void volume ratio (if any), and cure temperature. 

The following are the most important results from this investigation: 
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1. For unidirectional composites (ply), first-ply failure and laminate fracture were identical. Laminate theory was 
sufficient to evaluate behavior for uniaxial and combined loadings. 

2. For multi-ply laminates (angle-ply or multidirectional), laminate theory is sufficient to evaluate first-ply fail- 
ure. However progressive fracture is required to determine laminate fracture. 

3. In angle-ply and multidirectional laminates, residual stresses dramatically reduce first-ply failure in predomi- 
nant tensile stress states; however, their effects on laminate fracture stress are more benign. 

4. In comparison to uniaxial loading, combined tension and compression loadings had a negligible enhancement 
effect on laminate fracture stress. 

5. Angle-ply laminates loaded in biaxial tension exhibited substantial strain to fracture beyond the initial linear 
portion. This strain occurs with very small load increases. 

6. The pristine strength of the fibers was not exhibited in any of the laminates. Fiber fractures occurred at about 
70 to 80 percent of their respective pristine fracture strains. 
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TABLE L — MECHANICAL PROPERTIES OF VARIOUS MATRICES 


Property 

Matrix 


Hercules 3501-6 

DFVLR BSL914C 

Ciba Geigy 

Ciba Geigy 


epoxy 

epoxy 

LY556/HT907/DY063 

My750/HY917/DY063 






epoxy 

epoxy 


Pristine 

ICAN 

Pristine 

ICAN 

Pristine 

ICAN 

Pristine 

ICAN 

Modulus, GPa 

42 

5.8 

4.0 

5.724 

3.35 

4.772 

3.35 

4.2827 

Shear modulus, GPa 

1.56 

(a) 

1.481 

(a) 

1.24 

(a) 

1.24 

(a) 

Poisson’s ratio 

0.34 

0.38 

0.35 

0.398 

0.35 

0.4 

0.35 

0.4 

Strength, MPa 









Tensile 

69 

53.75 

75 

30.26 

80 

41.79 

80 

47.92 

Compressive 

250 

223.9 

150 

224.13 

120 

137.2 

120 

173.7 

Shear 

50 

94.068 

70 

93.79 

— 

86.89 


87.5 

Tensile failure strain, percent 

1.7 


4 

— 

5 


5 


Thermal coefficient, °C _1 

45x 10 6 

43.2 x 1 O' 6 

55x10* 

45x10* 

58x 10* 

64.8x10* 

58x10* 

52.38x10* 


a Shear modulus = (E/2(l + v)> is computed internally in ICAN code. 


TABLE II — MECHANICAL PROPERTIES OF VARIOUS FIBERS 


Property 

Fitx 

:r 


AS4 

T300 

H3HH9H 

Silenka E-glass 1200 tex 


Pristine 

ICAN 

Pristine 

ICAN 

Pristine 

ICAN 

Pristine 

ICAN 

Modulus, GPa 









Longitudinal 

225 

208.34 

230 

220.68 

80 

85.862 

74 

73.1 

Transverse 

15 

14.96 

15 

14.965 


83.860 

74 

81.72 

Shear modulus, GPa 









Inplane 

15 

14.96 

15 

14.965 

33.33 

33.31 


30.8 

Transverse 

7 

6.9 

7 

6.965 

33.33 

32.62 


30.8 

Major Poisson’s ratio 
Longitudinal strength, MPa 

02 

02 

02 

02 

02 

0.2 

02 

0.2 

Tensile 

3350 

3282.7 

2500 

2496 

2150 

1889 

2150 

2124 

Compressive 

Longitudinal failure strain, percent 
Tensile 

2500 

1.488 

6455 

2000 

1.086 

2034.48 

1450 

2.687 

1786 

1450 

2.905 

1310 





Compressive 

1.111 

0.868 

1.813 

1.959 





Thermal coefficient, °C“ I 





Longitudinal 

-0.5x10* 

-1.818x10* 

-0.7x10* 

-1.773x10-* 

4.9x10* 

6.633x10-* 

4.9x10* 

6.822x10* 

Transverse 

15x 10* 

16.182x10* 

12x10* 

12x10* 

4.9x10* 

4.86x10-* 

4.9x10* 

11.16x10* 


Property 


TABLE HI— MECHANICAL PROPERTIES OF VARIOUS PLIES 


Fiber/matrix 


AS4/Hereules T300/DFVLR 

3501-6 epoxy BSL914C epoxy 




Pristine 

ICAN 

Pristine 

ICAN 

Pristine 

ICAN 

Pristine 

ICAN 

Fiber content, vol% 
Modulus, GPa 

60 

60 

60 

60 

62 

62 

60 

60 

Longitudinal 

126 

127.309 

138 

140 

53.48 

53.37 

45.6 

45.634 

Transverse 

11 

11 

11 

11 

17.7 

17.57 

162 

16.206 

Inplane shear modulus, GPa 
Poisson’ ratio 

6.6 

6.289 

5.5 

6.248 

5.83 

6.427 

5.83 

5.937 

Major 

0.28 

0.272 

0.28 

027 

0.278 

0.28 

0.278 

0279 

Through thickness 
Longitudinal strength, MPa 

0.4 

0.414 

0.4 

0.44 

0.4 

0.403 

0.4 

0.415 

Tensile 

1950 

1969 

1500 

1500 

1140 

1131 

1280 

1295 

Compressive 

Transverse strength, MPa 

1480 

1480 

900 

900 

570 

570 

800 

800 

Tensile 

48 

48 

27 

27 

35 

35 

40 

40 

Compressive 

200 

200 

200 

200 

114 

114.48 

145 

144.964 

Inplane shear strength, MPa 
Thermal coefficient, °C _I 

79 

79.93 

80 

80 

72 

72.413 

73 

73 

Longitudinal 

-0.1x10* 

-0.1x10* 

-0.1x10* 

-0.940x10* 

8.6x10* 

8.6x10* 

8.6x10* 

8.64x10* 

Transverse 

26x 10* 

26x10* 

26x10* 

2624x10* 

26.4x10* 

26.26x10* 

26.4x10* 

26.28x10* 
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Figure 1 . — Simulation cycle of progressive fracture in composite laminates and structures via CODSTRAN. 



1 Equilibrium 

2 Incremental damage: 

properties degrade 

3 Additional damage: 

more degradation 

4 Equilibrium: no 

additional damage 

5 Next load increment 


Displacement 
Figure 2.— CODSTRAN load incrementation. 
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Load 



Figure 3. — Overall CODSTRAN simulation. 



Figure 4. — Upper half of ICAN-predicted biaxial failure 
envelope for 0° unidirectional E-Glass/LY556/HT907/ 
DY063 lamina, based on first-ply failure (laminate 
thickness = 1 mm). 



Figure 5. — ICAN-predicted biaxial failure stress enve- 
lope for 0° unidirectional T300/epoxy BSL914C 
lamina, based on first-ply failure (laminate thickness 
= 1 mm). 



a x , MPa 

Figure 6. — ICAN-predicted biaxial failure envelope 
of 0° unidirectional E-Glass/MY750/HY91 7/DY063 
lamina, based on first-ply failure (laminate thickness 
= 1 mm). 



ay, MPa 

Figure 7. — ICAN-predicted biaxial failure envelope 
of (90°/±30°) 2 E-Glass/HT907, based on first-ply 
failure (laminate thickness = 2 mm). 
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Figure 8. — Upper half of (CAN -predicted biaxial failure 
envelope of E-Glass/LY556/HT907, (907±30°) s , based 
on first-ply failure (laminate thickness = 2 mm). 



Figure 10. — ICAN-predicted biaxial failure envelope of 
(90°/+4570°) s AS4/Epoxy 3501-6, based on first-ply 
failure with and without residual stresses (laminate 
thickness = 1 .1 mm). 



oy, MPa 

Figure 9.-ICAN-predicted biaxial stress failure enve- 
lope of (±55°) s , E-Glass/M Y750/HY91 7/DY063, 
based on first-ply failure (laminate thickness = 1 mm). 



oy, MPa 

Figure 1 1 . — CODSTRAN-predicted failure envelope of 
0° unidirectional E-Glass/LY556/HT907/DY063 
lamina. 
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IE 


-300 -200 -100 


100 200 300 


oy, MPa 

Figure 12. — Stress in x-direction versus stress in 
/-direction for (907±30°/90°) E-Glass/LY556/HT907/ 
DY063 laminate (predicted by CODSTRAN code). 


« 150 



-300-150 0 150 300 450 600 750 900 

oy, MPa 

Figure 14. — Stress in x-direction versus stress in 
/-direction for ±55° E-Glass/M Y750/HY91 7/DY063 
angle-ply laminate (predicted by CODSTRAN code). 




c r x , MPa 

€y, percent 

Figure 1 3. — Shear stress versus transverse stress for 
(90°/±30790°) E-Glass/LY556/HT907/DY063 laminate 
(predicted by CODSTRAN code). 

Figure 15. — Stress-strain curve for E-Glass/M Y750/ 
HY917/DY063 ±55° angle-ply laminate (predicted 
by CODSTRAN code). 









, MPa -n <r y , MPa 



e x , €y, percent 


igure 16. — Stress in y-direction versus strains in 


x-direction for E-Glass/MY750/HY917/DY063 ±55° 


angle-ply laminate (predicted by CODSTRAN code). 



c x , Cy, percent 


Figure 18. — Stress versus x- and y-direction strains 
for AS4/3501-6 07±45790° quasi-isotropic laminate 
(predicted by CODSTRAN code). 



c y , percent 


Figure 17. — Stress versus y-direction strain for 
AS4/3501-6, 07±45790° quasi -isotropic laminate 
(predicted by CODSTRAN code). 


800 | — 



Cy, percent 


Figure 19. — Stress-strain in y-direction for E-Glass/ 
MY750/HY91 7/DY063 ±45° angle-ply laminate 
(predicted by CODSTRAN code). 
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e x ,€y, percent 


Figure 20. — Stress in y-direction versus strain in x- and 
/-directions for E-Glass/MY750/HY91 7/DY063 ±45° 
angle-ply laminate (predicted by CODSTRAN code). 



e x , percent 

Figure 21 . — Stress-strain curve for E-Glass/MY750/ 
HY917/DY063: 0°/90° cross-ply laminate (predicted 
by CODSTRAN code). 
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